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Abstract. We reconstructed the energy and the position of the shower maximum of air 
showers with energies E > 100 PeV applying a method using radio measurements performed 
with Tunka-Rex. An event-to-event comparison to air-Cherenkov measurements of the same 
air showers with the Tunka-133 photomultiplier array confirms that the radio reconstruction 
works reliably. The Tunka-Rex reconstruction methods and absolute scales have been tuned 
on CoREAS simulations and yield energy and X^ax values consistent with the Tunka-133 
measurements. The results of two independent measurement seasons agree within statistical 
uncertainties, which gives additional confidence in the radio reconstruction. The energy 
precision of Tunka-Rex is comparable to the Tunka-133 precision of 15 %, and exhibits a 20 % 
uncertainty on the absolute scale dominated by the amplitude calibration of the antennas. Eor 
Xmax) this is the first direct experimental correlation of radio measurements with a different, 
established method. At the moment, the X^ax resolution of Tunka-Rex is approximately 
40g/cm^. This resolution can probably be improved by deploying additional antennas and 
by further development of the reconstruction methods, since the present analysis does not yet 
reveal any principle limitations. 
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1 Introduction 

After more than 100 years of cosmic-ray measurements, the mass composition and the energy 
spectrum of the primary particles is relatively well-known only in the energy range accessible 
by direct measurements, E < 10^^ eV. At higher energies the flux of cosmic-rays is too low 
to perform direct measurements with statistical significance. Consequently, the knowledge 
is poorer and relies on the measurement of secondary-particle cascades called air showers. 
These indirect cosmic-ray measurements use the atmosphere as a calorimeter. The energy of 
the primary particle can be estimated from the energy contained in the shower. The type of 
the primary particle can, for example, be deduced from the longitudinal shower development. 
Heavy primary particles such as iron nuclei on average interact earlier in the atmosphere than 
light particles such as protons. Thus, the mass composition can be statistically estimated from 
measurements of the atmospheric depth of the shower maximum, Amax- 

Established air-shower techniques are the measurement of secondary particles on ground, 
the measurement of fluorescence or air-Cherenkov light by dedicated telescopes or non-imaging 
detectors on ground. The latter two methods have a relatively high accuracy for the energy 
and for Amaxj but are limited to dark and clear nights. Detection of the radio emission by air 
showers is an additional technique, which does not suffer from this intrinsic limitation of the 
exposure. The radio signal is sensitive to the shower energy [1] and A^ax [2“5], too, but it has 
not yet been demonstrated experimentally that the accuracy for both variables is competitive 
with those of the other techniques, as indicated by simulation-based studies, e.g., [6-8]. 
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Figure 1. Map of the antenna array Tunka-Rex and its host experiment, the photomultiplier array 
Tunka-133. 18 antennas had been available since the start of the first season in October 2012, and 
7 additional antennas have gone in operation only for the second season starting October 2013. The 
shower cores of the Tunka-Rex events are indicated for both seasons. The dashed circle denotes the 
cut on the central area with 500 m radius used for reconstruction. 


The main origin of the radio emission is the geomagnetic deflection of relativistic elec¬ 
trons and positrons in the shower, which induces a time-variable current [9]. The amplitude 
of this emission is proportional to the shower energy and to sin a, with a the angle between 
the geomagnetic field and the shower direction. On a weaker level the Askaryan effect con¬ 
tributes, i.e., radio emission due to the time variation of the net charge excess in the shower 
front [10-12]. The interference of both effects leads to an azimuthal asymmetry of the radio 
footprint [5, 13]. Finally, the refractive index of air causes Cherenkov-like effects for both 
emission mechanisms [14-16]. 

Tunka-Rex [17] is the radio extension of the Tunka observatory for cosmic-ray air showers. 
Its main detector, Tunka-133, is an array of non-imaging photomultipliers detecting the 
Cherenkov light emitted by air-showers in the atmosphere in the energy range of 10^® —10^® eV. 
Tunka-133 is fully efficient for all zenith angles Q < 50° at energies above 10^®'^ eV [18], which 
covers the full energy range of Tunka-Rex. For Tunka-133 measurements, the energy of the 
primary particle is reconstructed from the flux of the Cherenkov light at 200 m distance 
to the shower axis, and Xmax is reconstructed from the steepness of the amplitude-distance 
function. Using these methods, Tunka-133 features an energy resolution of about 15 %, and 
an Xjnax resolution of about 28g/cm^ [19]. For comparison, the average difference between 
the extreme cases of a pure proton and a pure iron composition is in the order of lOOg/cm^. 
Since Tunka-133 triggers Tunka-Rex, the same air showers are measured simultaneously with 
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Figure 2. Sky map of the Tunka-Rex events of both seasons. The efficiency rises with zenith angle 
and is smaller for arrival directions close to the geomagnetic field Bgeo- 


the air-Cherenkov and the radio detector. 

These combined measurements are used for a cross-check of both methods. In particular 
the precision of Tunka-Rex is estimated by a comparison of the energy and Xmax recon¬ 
structions to Tunka-133. Reconstruction methods for Tunka-Rex have been developed with 
CoREAS simulations [20] (the radio extension of CORSIKA [21]) taking into account previ¬ 
ous theoretical predictions [6, 7, 22, 23], and experience by other experiments, in particular 
LOPES [3], AERA [11] and LOFAR [24]. 

The parameters in the reconstruction methods have been determined from the simula¬ 
tions (cf. reference [8] for details), and not tuned against the Tunka-133 measurements. First 
results on energy and Xmax reconstruction have already been presented in reference [25] using 
the hrst season of Tunka-Rex measurements from October 2012 to April 2013 (effectively 280 
hours of measurements). The same reconstruction method is now applied also to the second 
season of data from October 2013 to April 2014 (effectively 260 hours of measurements), 
and the combined results and the cross-check between both measurement seasons are pre¬ 
sented in this article. In the experimental data we observe a direct correlation between the 
Tunka-Rex and Tunka-133 reconstructions, which is used to estimate precision and accuracy 
of Tunka-Rex for energy and Aimax- 

2 Experimental setup and data selection 

Tunka-Rex started in autumn 2012 with 19 antenna stations, and in summer 2013, 6 further 
antennas went into operation giving a total number of 25 antenna stations used for the present 
analysis. In the central area covering about Ikm^ the spacing is approximately 200 m. The 
outer antennas contribute only to a few events (see figure 1). The air-Cherenkov array Tunka- 
133 covers the inner area with 133 photomultipliers, and features additional photomultipliers 
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in the outer area. Upon a coincidence trigger of the photomultipliers, both the radio and 
the air-Cherenkov detector are read out in parallel. The electric field of the radio signal 
measured by Tunka-Rex is reconstructed in an effective bandwidth of 35 — 76 MHz. Details 
on the detector setup and its calibration can be read in reference [17]. 

The reconstruction of Tunka-133 is fully efficient and has a reliable reconstruction for 
events with zenith angles 9 < 50°. Therefore, only events with 9 < 50° have been selected for 
the present analysis (see figure 2). For these events the Tunka-133 reconstruction of the shower 
provides: energy of the primary particle atmospheric depth of the shower maximum 
^max) and axis, i.e., the direction and point of incidence on the ground (= shower core). For 
all Tunka-Rex events we used the shower core of Tunka-133 as input, because of the denser 
spacing of the photomultipliers compared to that of the antennas. The direction was only 
used as cross-check: events are considered false-positive, i.e., contaminated by background 
pulses and excluded from the analysis, if the directions reconstructed from the radio and 
air-Cherenkov arrival times deviate by more than 5° [17]. 

For the selection of Tunka-Rex events we applied a quality cut based on the signal-to- 
noise ratio in individual antenna stations. The signal has to exceed a signal-to-noise ratio, 
SNR = [S/N)"^ > 10 in at least three antenna stations, where the signal S is the maximum of 
an envelope on the electric-field strength of the pulse, and the noise N the root-mean square of 
the electric field in a time window before the signal. This cut is set such that pure background 
has a chance of about 5 % to pass in an individual antenna, which balances efficiency and 
purity. Given 25 antenna stations, this leads to a non-negligible probability that in a true 
event antenna stations with false-positive signals are contained. To exclude most of these 
false-positive signals, we sort the antennas by their distance to the shower axis. After two 
antenna stations failing the SNR cut, any further antenna stations are excluded from the 
analysis for the particular event. Thus, the risk that an event is contaminated by a station 
with false-positive signal is reduced to approximately 10 %. If surviving false-positive signals 
significantly impact the event reconstruction, they are removed by the cut requiring the 5° 
agreement between the direction reconstructed by Tunka-Rex and Tunka-133. Thus, the 
remaining data set is assumed to be practically free of false-positive signals. Finally, for the 
Aimax reconstruction we apply additional quality cuts: At least one antenna above threshold 
has to be at a distance Taxis > 200 m to the shower axis, since Aimax is reconstructed from the 
slope of the lateral distribution, which requires a sufficient lever arm. Moreover, the shower 
core has to be inside the central area of the array (dashed circle in figure 1), and the estimator 
for the statistical reconstruction uncertainty of Tunka-Rex must not exceed a certain value: 
CT(A-at) < 50g/cm2. 

Reconstruction of the electric-field strengths (= amplitude) in individual antenna sta¬ 
tions is performed with a modified version of the ‘Auger Offline’ software developed by the 
Pierre Auger Collaboration [26]. This software applies all known experimental characteristics 
to the recorded signals, in particular the frequency-dependent gain and pulse distortion of 
the various components in the signal chain as determined by calibration measurements [17]. 
For the reconstruction of the electric-field vector the simulated direction dependence of the 
antenna gain is used, which has been validated for a few selected directions by calibration 
measurements. 

For the suppression of narrow-band radio-frequency interferences (RFI) we apply rect¬ 
angular band-stop filters of ±0.2 MHz width (corresponding to 3 bins) at each integer 5 MHz 
(35, 40, ..., 75MHz), since these frequencies are often contaminated by RFI background 

^In earlier analyses we used a gliding median filter to suppress narrow-band RFI, but it turned out that this 
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Figure 3. Average ratio of the true pulse amplitude of CoREAS simulations ftrue and the amplitude 
after adding measured noise fmeas over the signal-to-noise ratio. 


Since the radio emission from the air shower is broadband, it is significantly less affected by 
the filter than the narrow-band disturbances. The background remaining after filtering still 
on average increases the measured pulse height Tmeas compared to the true pulse height Ttruej 
as already seen at the LOPES experiment [27] and assumed for early experiments [28]. We 
studied this effect by adding measured background to about 300 CoREAS simulations (fig¬ 
ure 3), and found that the parametrization given in reference [27] has to be slightly modified 
introducing a normalization factor A: = 4.1 adapted for Tunka-Rex. This normalization factor 
reflects the different ways how the noise level is measured in LOPES (mean of local maxima 
in the instantaneous amplitude) and Tunka-Rex (root mean square of the electric field, which 
corresponds to the mean power of noise). Thus, we use the following formula to correct for 
the average effect of noise: 


SNR’ 

with the signal-to-noise ratio SNR defined above. 

An overview of the available event statistics is provided in table 1, and figure 4 shows 
the distribution of the selected events over energy and X^ax- Most of the events used for the 
energy reconstruction are at energies above 10^^ eV. The fraction does not reach 100 % because 
some Tunka-133 events have small geomagnetic angles. Also, there were occasional technical 
issues in single antenna stations. For the high-quality events used for the X^ax reconstruction, 
the threshold is slightly higher: above an energy of about 6 • 10^^ eV, all events used for energy 
reconstruction pass also the high-quality cuts required for X^ax reconstruction (see figure 5). 
Apart from technical issues (e.g., missing or malfunctioning antennas), the fraction shown 
in figure 5 is the efficiency of Tunka-Rex for zenith angles 9 < 50°, because the reference 
is the fully efficient Tunka-133 trigger. The total number of Tunka-Rex events passing the 
quality cuts is larger [17], since in contrast to the air-Cherenkov detection, radio detection is 

filter significantly increases the impact of background on the measured amplitudes, without significantly in¬ 
creasing the signal-to-noise ratio. Nevertheless, the exact method of RFI suppression has only minor influence 
on the high-quality events used in the present analysis. 
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Table 1. Statistics of Tunka-Rex events triggered at zenith angles 9 < 50° in both seasons after 
certain quality cuts. The cuts are applied consecutively in the given order. The event numbers denote 
the number of surviving events after application of all cuts up to that point. The first three cuts are 
used to reject false-positive events contaminated by background pulses. The remaining 91 -I- 87 events 
are used for energy reconstruction. For reconstruction three additional quality cuts are applied 

leaving 23 -I- 19 events for the two measurement seasons, respectively. 


Number of events 
after cut on 

first season of 

Oct. 2012 - Apr. 2013 
(280 effective hours) 

second season of 

Oct. 2013 - Apr. 2014 
(260 effective hours) 

> 3 stations with signal 

244 

445 

rejecting outlier stations 

122 

147 

Z (Tunka-Rex, Tunka-133) < 5° 

91 

87 

at least one antenna at Taxis > 200 m 

64 

56 

a(A-t) < 50g/cm2 

25 

22 

inner area (r < 500 m) 

23 

19 




reconstructed atmospheric depth of shower maximum (g/cm*) 


Figure 4. Distributions of energy, Flpr, and atmospheric depth of the shower maximum, X„iax, of the 
Tunka-Rex events used in the present analysis. The figure shows the reconstructed values by both 
Tunka-Rex and Tunka-133. Bars represent the statistical uncertainty calculated as \fN. 


more efficient for inclined air showers. Nevertheless, for the present analysis only events with 
6 < 50° have been selected, since Tunka-133 is fully efficient for this zenith angle range. 

3 Reconstruction of energy and shower maximum 

In the present analysis, the reconstruction of energy and shower maximum is based on the 
lateral distribution, i.e., the dependence of the radio amplitude on the distance to the shower 
axis. This approach has been used earlier. Different experiments have shown that the am¬ 
plitude at a certain distance is correlated with the primary energy, after correcting for the 
geomagnetic angle [1, 29-31]. Moreover, it has been shown that the shape of the lateral 
distribution depends on the distance to the shower maximum [2, 3, 5]. 

Our reconstruction method in brief is: 

1. We correct the lateral distribution for the azimuthal asymmetry caused by the interfer¬ 
ence of the geomagnetic and the Askaryan effects, and for the geomagnetic angle a. 
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Figure 5. Fraction of events which passes the quality cuts applied for the Tunka-Rex energy and Amax 
reconstructions, respectively, in comparison with all events triggered by the fully efficient Tunka-133 
array in the inner area (r < 500 m). 


2. We fit a Gaussian LDF, the shape of which depends on three parameters: One scale 
parameter, S 120 , which is the amplitude at 120 m, and two shape parameters ai and 02 , 
which describe the slope and the width of the Gaussian LDF, respectively. £120 and ai 
are ht to the data points, 02 is fixed using a parametrization depending on energy and 
zenith angle. 

3. We use the fit result £120 to reconstruct the primary energy and the fit result for ai 
at 180 m axis distance, representing the slope of the lateral distribution, to reconstruct 
the atmospheric depth of the shower maximum X^ax- 

3.1 Asymmetry correction 

Because of the different polarizations of the radio emission generated by the Askaryan and 
the geomagnetic effects [11, 12, 32], the radio amplitude on ground depends not only on the 
geomagnetic angle a and distance to the shower axis r, but also on the azimuth angle (jig 
of each antenna relative to the shower axis. Fitting the asymmetry for individual events by 
a two-dimensional lateral distribution function (LDF) introduces additional parameters [5]. 
To maintain events with low station multiplicity, we have developed an alternative approach 
described in reference [8], which is well-suited for sparse arrays like Tunka-Rex. 

Our approach makes use of the fact that the asymmetry to first order only depends 
on the shower geometry, which can be reconstructed independently from the LDF by arrival 
time measurements. Based on GoREAS simulations made for the Tunka-Rex setup, the 
azimuthal asymmetry of the footprint is parametrized as function of the antenna position 
using a constant typical value of Easym = 8.5 % for the relative strength of the Askaryan 
effect.^ Using this parametrization, we correct all measured amplitudes in individual antennas 

^ Although the size of the asymmetry Easym depends on zenith angle and distance to the shower axis [12], 
simulation studies for Tunka-Rex have shown that parameterizing these dependences of Easym instead of 
using a constant value does not describe the lateral distribution better. Probably this is because Tunka-Rex 
measurements mostly fall in a limited zenith range. Moreover, Easym shows little variation in the distance 
range of 100 — 200 m used for energy and Amax reconstrnction. 
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Figure 6. Example event. Left: footprint. Large crosses are triggered stations, and blue crosses are 
those stations with signal above the signal-to-noise threshold. The small dark-gray cross is a station 
without measurement for this event. The star and the line indicate the reconstructed shower core 
and direction. Right: Lateral distribution before (open circles) and after (filled circles) asymmetry 
correction, but with correction for the geomagnetic angle in both cases. The LDF is fitted to the 
corrected points. The light gray point is a station with signal below the signal-to-noise threshold. 


for the azimuthal asymmetry, and obtain a one-dimensional lateral distribution in which the 
amplitude depends only on the distance to the shower axis. Simultaneously, we correct for 
the relative strength of the geomagnetic Lorentz force, since this depends only on the shower 
geometry, too, namely on the geomagnetic angle a. Thus, we use the following equation 
to obtain the corrected amplitude Scovif), which is later used to fit a radially symmetric 
lateral-distribution function: 


£cor{r) = Sir, (j)g)/^J+ 2easym cos (pg sin a + sin^ a , (3.1) 

with £{r,(j)g) the measured amplitude after noise correction (£^true of equation 2.1), r the 
distance to the shower axis, cjig the azimuth angle relative to the shower axis (with (/>g = 0, 
when a station is in the direction of the geomagnetic Lorentz force), a the geomagnetic angle, 
Easym = 0.085 the size of the asymmetry, and £’cor(^) the amplitude after correction for the 
asymmetry and the geomagnetic angle. After this correction of the amplitude measured at 
each antenna station, the lateral distribution is approximately azimuthally symmetric around 
the shower axis. Thus, to first order the amplitude depends only on the distance to the 
shower axis r, and the lateral distribution can be fitted with a one-dimensional function the 
parameters of which depend on energy and distance to the shower maximum. 

Figure 6 shows an example event with many stations, where the asymmetry correction 
leaves the energy reconstruction almost unchanged, since the LDF fit averages the asymmetry 
out. However, the Xmax reconstruction is slightly improved for this event (Tunka-Rex is 

613±15g/cm^ and 621±15g/cm^ before and after correction, respectively, and the Tunka-133 
Xma.x is 657 ± 28g/cm^). Especially for events with few stations, the asymmetry correction 
is important, and can have an effect on the energy reconstruction in the order of 10 %. 










3.2 Lateral distribution function (LDF) 

The amplitude (= electric field strength) of this one-dimensional lateral distribution is propor¬ 
tional to the shower energy, because the radio emission in the shower is coherent. Moreover, 
the shape of the lateral distribution depends mainly on the distance to the shower maximum 
and has two main features: first, a bump caused by Cherenkov-like effects [14, 15, 32, 33], 
the existence, size and position of which depend on the distance to the shower maximum; 
second, an exponential tail [1, 34], the slope of which depends on the distance to the shower 
maximum, too. This means that even a one-dimensional LDF should contain at least three 
free parameters: one for the energy scale, and two for the shape, namely the width of the 
bump, and the slope of the exponential fall-off. For this purpose we use the following lateral 
distribution function: 


Scor{r) = Sro exp (ai(r - ro) + a 2 {r - r^f) . (3.2) 

with the amplitude Tcor(^) after correction for the azimuthal asymmetry and for the geomag¬ 
netic angle as function of axis distance r, and the fitted scale parameter which is the 
amplitude at the axis distance rg. The shape parameter oi determines the slope of the expo¬ 
nential fall-off, and the shape parameter 02 is related to the width of the Cherenkov bump. 
The distance parameter rg is a purely technical parameter. Using a minimizing chi-square 
fit, the shape of the function does not depend on the choice of rg, i.e, for different rg the 
resulting function is exactly the same, but just described with a different set of values for the 
parameters <^1 and 02 . However, the choice of rg affects the uncertainties and correlations 
of the fit parameters, which is why a specific choice of rg can simplify the reconstruction of 
Fpr, Xmax) and their statistical uncertainties. We set this distance parameter rg to 120 m 
for the reconstruction of the primary energy based on T 12 O) and to 180 m for the X^ax 
reconstruction based on oi, as these turned out to be the optimum values in the CoREAS 
simulation study. 

In the CoREAS simulations made for Tunka-Rex we observe that the scale parameter 
T 120 is strongly correlated with the primary energy Fpr. The shape parameter ai, i.e., the slope 
at 180 m axis distance is strongly correlated with the distance to the shower maximum. The 
second shape parameter 02 is correlated with ai, and weakly depends on energy and zenith 
angle. Thus, fitting 02 as a free parameter would complicate the reconstruction of Xmax 
significantly. For this reason we decided to perform the fit with only two free parameters, 
and to fix the parameter 02 using a parametrization depending on zenith angle 9 and the 
primary energy Fpr (details in appendix A). After fixing 02 , the LDF is fit to the data 
points. Ti 2 g carries the whole sensitivity to the primary energy Fpr, and the remaining free 
shape parameter oi carries the whole Xmax sensitivity of the LDF. Furthermore, reducing 
the number of free parameters in the LDF from three to two brings the additional advantage 
that we can use events with only three antenna stations above threshold (about half of the 
Tunka-Rex events in the present analysis). 

3.3 Primary energy Fp^, and atmospheric depth of shower maximum A m ax 

Energy Fpr and Amax are determined using the equations presented in reference [8]. While 
for Amax we use exactly the same equation as in the reference, i.e., the same values for all 
parameters in the equation, for Fpr we use a simplification. The simplification is to assume 
that the energy is exactly proportional to the amplitude, as expected for 100% coherent 
emission. This means that the exponent b = 0.93 of reference [8] is set to 1 instead and, thus, 
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Figure 7. Correlation and relative deviation (= difference divided by average) of the shower energy 
reconstructed with Tunka-Rex radio and Tunka-133 air-Cherenkov measurements. 


is not present in the equation used here. The reason for applying this simplification is that with 
the present statistics and quality of the measured data, we have found no evidence for b ^ 1. 
Moreover, there is no significant difference in the resulting energy precision whether b is set 
to 0.93 or 1. Consequently, at the moment we lack experimental justification for introducing 
b as an additional parameter and decided for the simpler solution. The simplification requires 
a different proportionality coefficient kl, which has been determined with the same set of 
CoREAS simulations. Summarizing, we use the following equations for the reconstruction of 
Epi- and Aimax with Tunka-Rex: 


Epr = Kl ■ £i20 , (3.3) 

-^max = X^et/cos9 - {A +Blog{ai+ b)), (3.4) 

where = 955g/cm^ is the atmospheric depth of the detector (typical value for the 

site used in the simulations), and the parameters are determined with CoREAS simulations: 
Kl = 884EeV/(V/m), A = —1864 g/cm^, B = —566g/cm^, oi given in units of meters, and 
b = 0.005. As explained above, the technical parameter of the LDF, lq, is set to 120 m for 
the determination of the amplitude at 120m, £ 120 , and to 180m for the determination of the 
slope parameter ai. 

The energy and X^ax reconstructions based on simulations were already presented in ref¬ 
erence [8]. The following section presents the experimental result of Tunka-Rex measurements 
and the comparison to the Tunka-133 air-Cherenkov measurements. 

4 Results 

There is a clear correlation between the energy reconstructed from the radio amplitude mea¬ 
sured by Tunka-Rex and the energy reconstructed from the air-Cherenkov light measured by 
Tunka-133 (figure 7). This indicates that the energy reconstruction works reliably with the 
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Figure 8. Correlation and difference of the distance to Xmax (atmospheric depth between the shower 
maximum and the shower core) as reconstructed with Tunka-Rex radio and Tunka-133 air-Cherenkov 
measurements. 

presented method. When combining both seasons, the average deviation between the Tunka- 
Rex and Tunka-133 energy values is (17 ± 2) %, which is only slightly larger than the 15% 
energy resolution reported by Tunka-133 [19]. 

Also for the position of the shower maximum there is a clear correlation between the 
radio and air-Cherenkov reconstructions (figure 8). Due to the stricter quality cuts the event 
statistics is lower than for the energy correlation. Still, the correlation is significant. This 
again indicates that also the reconstruction of Xmax works reliably for the selected, high- 
quality radio events. Because of the lower statistics, the A^max reconstruction precision of 
Tunka-Rex cannot be estimated very accurately. The standard deviation of the A"max differ¬ 
ence between both detectors is (47 ±5) g/cm^, which reflects the combination of the unknown 
Tunka-Rex uncertainty and the Tunka-133 uncertainty of 28g/cm^ [19]. This corresponds 
to a radio-only precision of better than 40g/cm^, when assuming that the resolutions of 
Tunka-133 and Tunka-Rex add in squares to the total deviation. 

For both the energy and Xmax reconstructions the observed correlation is statistically 
compatible with a 1:1 correlation of Tunka-133 and Tunka-Rex: the means of the Gaussians 
fitted in figures 7 and 8 are compatible with 0, where 0 corresponds to no bias (= systematic 
scale offset). This means that the energy and Xmax scales of both detectors agree on an 
absolute level within the measurement uncertainty. This is remarkable, since the Tunka- 
Rex scale is defined by CoREAS simulations and has not been tuned against Tunka-133. 
With already consistent scales there is no need for a true cross-calibration, i.e., en explicit 
calibration of both detectors against each other. 

5 Discussion 

5.1 Validity of the results 

The present analysis is the first event-by-event comparison of radio and air-Cherenkov mea¬ 
surements of the same air showers. Since both detection methods are sensitive to the energy 


11 



















and the position of the maximum of the electromagnetic shower component, a strong correla¬ 
tion is expected between radio and air-Cherenkov observables. This is especially noteworthy 
since the Tunka-Rex reconstruction is based on simulations and not tuned against the Tunka- 
133 measurements we compared with. Only details of the reconstruction method have been 
decided after looking at the data of the first season, e.g., the filter used to remove narrow-band 
interferences or the exact quality cuts. Neither of these small decisions made the and 

energy correlations appear or vanish. 

As additional cross-check of the validity of the results, we decided to blind the energy 
and Xmax reconstruction of Tunka-133 for the second season. These values were revealed 
to us only after we froze the reconstruction method, including digital filters and quality 
cuts, and after we decided to use the simplified equation 3.4 for energy reconstruction. In 
appendix B the results of both seasons are compared, and found to be compatible within 
statistical uncertainties. Thus, we conclude that our results are not corrupted by any kind 
of implicit tuning on the data set of the first year. Moreover, the quality cuts are general 
enough to cover the slightly different array configurations in both seasons, i.e., the missing 
antennas in the first season do not have significant impact on any of the present results. 

Potential biases not excluded by the cross-check concern possible special events. For 
example, events with their shower maximum very close to the detector will have a very steep 
lateral distribution and a small footprint. Thus, they will likely not be detected or will be 
rejected by the quality cuts. Still, since Tunka-133 is fully efficient and has a reliable recon¬ 
struction in the full zenith and energy range used in this study, the non-detection of certain 
events is not assumed to cause any significant bias in the observed correlations. Furthermore, 
at high energies less than a third of the events are rejected (cf. figure 5), and this due to 
technical reasons (missing antenna stations) or due to small geomagnetic angles leading to 
small radio amplitudes. Thus, the potential problem of a bias introduced by the rejection of 
any kind of special event can concern only a small fraction of events. Nevertheless, this has 
to be studied in more detail when larger statistics will be available. 

Consequently, the cross-check of both seasons cannot totally exclude all imaginable 
systematic biases implicit in the experimental setup or in the simulations simulations, but it 
definitely confirms that the observed correlations are real. 

5.2 Precision and accuracy 

The spread of the deviation between the Tunka-Rex and Tunka-133 reconstruction values for 
the same events can be used to estimate the reconstruction precision of Tunka-Rex. Assuming 
independent Gaussian uncertainties of both reconstructions, the standard deviation of the 
difference or ratio between Tunka-Rex and Tunka-133 is the quadratic sum of the absolute or 
relative uncertainties, receptively. If Tunka-Rex would feature the same precision as Tunka- 
133, thus, this standard deviation should be a factor of \/2 larger than the Tunka-133 precision 
alone, i.e.: \/2-15 % = 21 % for the energy, which is slightly larger than the observed standard 
deviation of (17± 2) %. Consequently, the energy precision of Tunka-Rex seems to be at least 
equal to that of Tunka-133. 

For Amax) we would expect a standard deviation of the difference of \/2 • 28g/cm^ = 
40g/cm^ if Tunka-Rex had the same precision as Tunka-133 (28g/cm^), but observe a stan¬ 
dard deviation of the difference of (47±5) g/cw?. The derived Tunka-Rex precision is slightly 
better than 40g/cm^, as expected from CoREAS simulations when including realistic back¬ 
ground [8]. Since the average spacing of the Tunka-133 photomultipliers is half of the average 
antenna spacing, this result does not necessarily imply that the radio Amax precision is in- 
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trinsically worse than the air-Cherenkov precision. Instead, that later studies will show how 
the Tunka-133 and Tunka-Rex resolutions compare for equal detector spacing, and by how 
much the Tunka-Rex precision will improve due to additionally deployed antennas. 

Using CoREAS simulations we searched for systematic uncertainties potentially degrad¬ 
ing the resolution further, in particular biases of the Xmax reconstruction over energy or zenith 
angle, but we did not find any significant effect. The hnally achievable Amax uncertainty of 
the radio technique might become equal to the one of the air-Cherenkov technique, when 
events of even higher quality are used. The spread between the Tunka-133 and Tunka-Rex 
Amax values approximately corresponds to the cut placed on the estimator for the 

statistical uncertainty of Tunka-Rex. This means that this estimator is reliable and can be 
used to set stronger quality cuts when higher statistics are available. 

In addition to the precision also the total accuracy is of importance. This is dominated 
by the 18% scale uncertainty of the amplitude calibration [17], but at least three further 
systematic scale uncertainties have to be considered: 

First, if the fraction of the total primary energy going in the electromagnetic shower 
component was simulated wrongly in CORSIKA, the energy scale would be wrong by roughly 
the same value. However, the electromagnetic shower component is relatively well under¬ 
stood, and as estimated by other experiments [35], this uncertainty on the ‘invisible’ energy 
is only a few percent. Moreover, we assume that the CoREAS radio extension of CORSIKA 
simulates the absolute radio amplitude emitted by the electromagnetic shower component 
correctly, i.e., with errors which are small against the calibration scale uncertainty of 18%, 
which is supported by recent experimental tests [17, 36]. Still, the fraction of energy go¬ 
ing in the electromagnetic shower component depends on the mass of the primary particle. 
Thus, the parameters in the method slightly depend on the assumed mass composition of 
the primary particles, in particular also the parameter kl determining the energy scale. For 
the present analysis, we determined kl as average of proton and iron simulations with the 
hadronic interaction model QGSJET-11.04 [37]. If we had taken a pure proton or a pure iron 
composition to determine k,l of equation 3.4, the reconstructed energy would be 4% smaller 
or larger, respectively. 

Second, there is some freedom in the equation used for energy reconstruction. For 
example, if we had used exactly the energy equation of reference [8] instead of equation 3.4, 
then the reconstructed primary energy would be about 10 % lower, although both equations 
have been tuned against the same set of CoREAS simulations. This difference is well within 
the 18% scale uncertainty of the amplitude calibration, but indicates that a more detailed 
study of the reconstruction method will become necessary once a more accurate calibration 
becomes available. 

Third, using CoREAS simulations we found a systematic bias of the energy reconstruc¬ 
tion depending on the distance to the shower maximum. The mean effect is smaller than 
the difference between proton- and iron-initiated showers, and thus negligible. However, for 
very close or far distances to the shower maximum (< 350g/cm^ or > 900g/cm^) the bias 
becomes larger than 10%, and will have to be corrected for in such rare events. Conse¬ 
quently, all known systematic uncertainties are small compared to the scale uncertainty of 
the amplitude calibration of 18 %. Thus, the total accuracy of the Tunka-Rex energy scale is 
approximately 20%. 

The Tunka-Rex energy scale defined by CoREAS simulations agrees within the un¬ 
certainties to the Tunka-133 scale. This means that the measurements of Tunka-Rex and 
Tunka-133 are not only correlated, but in addition agree on an absolute level, which is an- 
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other confirmation, that CoREAS seems to predict the radio signal correctly. 

5.3 Comparison with other experiments 

The 15 % energy precision of Tunka-Rex is at least as good as that of other radio arrays, like 
AERA [38], LOPES [3], and CODALEMA [39], which all reported energy precisions around 
20%. The Tunka-Rex scale uncertainty of 20% roughly corresponds to that of AERA [38] 
and LOPES [36], which compared the radio measurements to air-fluorescence and particle 
measurements on ground. LOFAR as well features an absolute calibration of the radio ampli¬ 
tude with similar accuracy [40], but experimental checks of the energy accuracy are limited 
by the poor energy resolution of the LORA particle detector array [5]. For X^ax reconstruc¬ 
tion, the resolution of LOPES has been significantly worse, using a simpler method in a more 
radio-loud environment than Tunka-Rex. The much denser array LOFAR is located in a more 
radio-quiet environment similar to the situation of Tunka-Rex. Its X^ax reconstruction is also 
based on CoREAS simulations using two different methods. With an X^ax reconstruction 
based on a htted LDF (like in our approach), the precision of X^ax at LOFAR is only slightly 
better than our resolution [5]. However, with a more computing-intensive method using many 
simulations for each individual event [4], the precision by LOFAR is at least twice as good 
as ours. This indicates that in addition to the planned deployment of additional antennas, 
also further improvements in the reconstruction methods could increase the Xmax resolution 
of Tunka-Rex. 

The most important scientific value of the present analysis is the experimental compari- 
sion of radio measurements against another, established technique. For the first time, energy 
and Xmax reconstructions based on absolutely-calibrated radio measurements have been com¬ 
pared to air-Cherenkov measurements. While the principle sensitivity of the radio signal on 
the longitudinal shower development was already demonstrated experimentally by comparing 
LOPES lateral distributions to measurements of the KASCADE muon-tracking detector [2], 
KASCADE did not feature direct X^ax measurements. Thus, with the present comparison 
of Tunka-Rex and Tunka-133, for the first time the radio reconstruction of X^ax is directly 
cross-checked with an independent X^ax measurement. This also gives more confidence in 
the results of other radio arrays whose reconstruction procedures are developed with the same 
CORSIKA + CoREAS Monte Carlo codes. 

6 Outlook 

The presently achieved accuracy for energy and X^ax reconstruction by Tunka-Rex is not yet 
at its principle limits. In particular for the X^ax reconstruction, there is room for improve¬ 
ment. As shown by LOFAR [4], the precision can be as good as 20 g/cm^, when using a denser 
radio array and a simulation-driven method. Thus, we plan to deploy additional antennas 
and to further improve the reconstruction method. With an increased antenna density, the 
average event will contain more antenna stations with signal, which allows fitting of lateral 
distributions with more free parameters. For example, the weak energy and X^ax sensitivity 
of the LDF parameter a 2 could be exploited when more antenna stations contribute to the 
fit. 

Compared to this study the antenna density will be tripled. At each of the 19 Tunka 
clusters in the inner area, already now a second antenna station is deployed, and a third one 
will be deployed in summer 2016. This enables an experimental test how the antenna density 
and the core resolution affect the X^ax resolution. Moreover, additional methods for X^ax 
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reconstruction based on other quantities of the radio signal can be used to improve the total 
accuracy. The shape of the wavefront [41] can be reconstructed by arrival time measurements, 
and the slope of the frequency spectrum [42] by sampling the pulse shape. Future studies 
will show by how much these methods can increase the total accuracy under practical 

conditions. 

The statistics of Tunka-Rex will be dramatically increased by joint measurements with 
the newly deployed particle-detector array Tunka-Grande [43], which provides a day-time 
trigger for Tunka-Rex. Tunka-Grande consists of 19 stations, one at each inner cluster, where 
each station features scintillator detectors on ground and under ground, for measurement of 
secondary air-shower electrons and muons, respectively. While unimportant for the present 
analysis focused on the feasibility of the reconstruction methods, for reconstruction of the 
cosmic-ray composition possible systematic uncertainties like selection biases are relevant. Us¬ 
ing the additional statistics and independent measurements provided by Tunka-Grande such 
systematic uncertainties can be checked experimentally. Moreover, the additional muon mea¬ 
surements can enhance the total accuracy for the mass composition, since the electron-muon 
ratio provides complementary mass information to Xmax [44]. In addition, the Tunka-Grande 
measurement accuracy will be enhanced by a cross-calibration on hybrid measurements with 
Tunka-Rex. 

Finally, the absolute energy measurement of air-showers with a radio array opens prospects 
to cross-calibrate the energy scales of different air-shower experiments via radio measurements. 
This aspect is worth to be investigated more deeply, in particular since antenna arrays can 
be a very economic add-on to existing detectors as shown by Tunka-Rex. 

7 Conclusion 

We have compared energy and X^ax reconstructed from Tunka-Rex and Tunka-133 mea¬ 
surements of the same air showers in the energy range Fipr > 10^^ eV. The reconstruction 
methods have been developed and tuned using simulations with GORSIKA and its radio ex¬ 
tension GoREAS. For both parameters we find a strong correlation between Tunka-Rex and 
Tunka-133 consistently in two seasons of data taking. This confirms experimentally that the 
radio measurements are indeed sensitive to energy and shower maximum. For X^ax this is 
the first direct confirmation based on a cross-check with a different, established experimental 
technique, namely air-Gherenkov measurements. 

The Tunka-Rex energy precision seems to be at least as good as the published Tunka-133 
resolution of 15 %. The total scale uncertainty of Tunka-Rex is dominated by the uncertainty 
of the amplitude calibration, and in total is in the order of 20 %. This is comparable to 
the scale uncertainty of particle detector arrays, like KASGADE-Grande [45]. Since the 
calibration uncertainty dominates the total scale uncertainty of Tunka-Rex, any improvement 
in the calibration accuracy will directly propagate to the energy-scale accuracy. Gonsequently, 
further efforts on the calibration will be necessary to increase the accuracy to the same level 
as the currently leading fluorescence and air-Cherenkov techniques. At the Pierre Auger 
Observatory, the scale accuracy of fluorescence measurements is 14% [35]. This now is in 
reasonable reach for future radio measurements, which unlike fluorescence measurements are 
available around-the-clock. Hence, radio measurements can be used to determine the absolute 
energy scale of air-shower measurements. 

The Xmax precision of Tunka-Rex is roughly 40g/cm^, and can be slightly increased by 
setting stricter quality cuts at the cost of statistics. This resolution is sufficient to statistically 
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distinguish light from heavy primary particles. The resolution is worse than the one currently 
achieved by non-imaging air-Cherenkov measurements (28g/cm^ for Tunka-133 [19]) and by 
air-fluorescence measurements (20 g/civ? for the Pierre Auger Observatory [35]). Nevertheless, 
the resolution of Tunka-Rex is not yet at its limit and can likely be increased by further 
improvements of the reconstruction method, and by deploying additional antennas. 

For astrophysical applications the total accuracy for the reconstruction of the mass 
composition as a function of energy counts, e.g., to better study the transition from Galactic to 
yet-unknown extra-galactic cosmic-ray sources assumed in the energy range of Tunka-Rex [46]. 
Tunka-Rex will provide additional statistics exactly in the energy range where current Tunka- 
133 analyses are limited by statistics. Moreover, hybrid measurements of air-showers by 
Tunka-Rex and Tunka-Grande started in autumn 2015. These hybrid measurements can be 
used to enhance the total accuracy for the mass composition by combining radio and muon 
measurements. 
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A Parametrization of LDF parameter 02 

In the Gaussian lateral distribution function (LDF, equation 3.2), the parameter 02 defining 
the bump has been set by the following parameterization determined with CoREAS simula¬ 
tions made for the situation of Tunka-Rex [8]: 
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Since the 02 parameter is a function of the primary energy Fpr, the shape of the Gaussian 
LDE used to reconstruct the primary energy depends implicitly on the energy itself. This 
reciprocal dependence on can either be solved by an iterative approach, or by using a 
pre-estimate for Ep^-. We decided for the latter solution, using the energy estimator based 
on a simpler exponential LDF presented in reference [17], since its precision is only slightly 
worse than than in our more advanced approach here. 


B Comparison of both seasons 

The data set of Tunka-Rex is split in two seasons of about equal size. As additional cross¬ 
check of the validity of our methods, the Tunka-133 reconstruction of energy and Vmax had 
been blinded for the second season. Only for the first season the Tunka-133 energy and Vmax 
values had been known, and were used for occasional cross-checks during the development 
of the Tunka-Rex methods. After we frozen the Tunka-Rex reconstruction methods, we 
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Figure 9. Comparison of the first and second season: reconstructed energy (cf. figure 7). 


predicted the energy and for the second season based on the radio measurements. Only 

afterwards, Tunka-133 unblinded their reconstruction of the second season to us. 

The plots in this appendix show the comparison of the Tunka-Rex and Tunka-133 energy 
and Xmax reconstructions (hgures 9 and 10) separately for both seasons. The results of both 
seasons are consistent within statistical uncertainties. The quality cut on events inside of the 
inner area (r < 500 m, dashed circle in hgure 1) was discovered a posteriori by us, i.e., after 
we had frozen the reconstruction methods for the unblinding procedure. Thus, in figure 10 5 
additional events are present (2 in the first season, 3 in the second season), which are outside 
of the inner area and slightly degrade the X^ax resolution. 

As supplementary material we will upload the list of the events used for the present 
analysis. This list contains the reconstructed energy and X^ax values of the events. However, 
the values should not be used for reconstruction of the cosmic-ray energy spectrum or the 
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Figure 10. Comparison of the first and second season: distance to shower maximum (cf. figure 8). 


mass composition. While unimportant for the present study comparing the Tunka-Rex to the 
Tunka-133 reconstruction, for such analyses selection biases have to be taken into account, 
which cannot be derived from the event list alone. 

The original event lists for the first ‘tuning’ and second ‘prediction’ seasons, as created 
after freezing the reconstruction methods, but before unblinding, are available at: 
http://www.ikp.kit.edu/tunka-rex/ 

They can be decrypted using the following Linux command: 
openssl aes-256-cbc -d -in encryptedFile -out decryptedFile 
The passwords are ‘TunkaRexl_H9AoxFywAt’ for the tuning season, 
and ‘TunkaRex2_z3DumFNfsq’ for the prediction season. 

All these events are also contained in the above-mentioned list submitted as supplementary 
material. 
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